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Abstract 

We present an approach to derive longitudinal shower development information from the longitudinally unsegmented 
BABAR electromagnetic calorimeter by using tracking information. Our algorithm takes advantage of the good three- 
dimensional tracking resolution of BABAR, which provides an independent geometric constraint on the shower as 
measured in the BABAR crystal calorimeter. We show that adding the derived longitudinal shower development 
information to standard particle identification algorithms significantly improves the low-momentum separation of 
pions from electrons and muons. We also verify that the energy dependence of the electromagnetic shower development 
we measure is consistent with the prediction of a standard electromagnetic shower model. 
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1. Introduction 

In present-day nuclear and particle physics ex- 
periments, inorganic scintillating crystals, such as 
Nal(Tl) and CsI(Tl), are often used to construct 
electromagnetic calorimeters when a precise mea- 
surement of the energy is required [1]. Crystal 
calorimeters can be finely segmented in the dimen- 
sion transverse to the shower development without 
sacrificing energy resolution, thus providing a good 
measurement of the lateral shower development. 
However, engineering and energy resolution consid- 
erations prevent finely segmenting crystal calorime- 
ters along the direction of shower development. 
Both lateral and longitudinal shower development 
information are useful in charged particle identi- 
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fication (PID) algorithms, particularly in electron 
identification. Because crystal calorimeters cannot 
provide direct longitudinal shower development in- 
formation, they lose an important input to particle 
identification. 

In this paper, we present a technique in which 
longitudinal shower development information is 
indirectly extracted from a longitudinally unseg- 
mented crystal calorimeter in conjunction with a 
precise tracking system. This technique was devel- 
oped for use with the BABAR. [2] detector, but it can 
be applied at any detector which combines crystal 
calorimetry and precision tracking. It exploits the 
fact that BABAR has a tracking system capable of 
precisely determining the three-dimensional trajec- 
tory of charged particles, and the fact that these 
trajectories are not in general collincar with the 
crystal axes. A similar algorithm has previously 
been used for sampling calorimeters with fine lateral 
segmentation [3]. 
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We also explore the usefulness of the derived lon- 
gitudinal shower information for particle identifi- 
cation at BABAR. We show that using this infor- 
mation as part of an electron identification algo- 
rithm improves the electron vs. pion discrimination, 
particularly at momenta below 600 MeV/ c. We also 
show that this information can be generally used 
to improve the separation of the five most common 
charged particles (e, /i, it, K and p) , owing to the dif- 
ferent longitudinal shower development in a crystal 
calorimeter of these different particle types. In par- 
ticular, we find that the longitudinal shower infor- 
mation significantly increases the minimum momen- 
tum reach at which muons can be separated from 
pions in BABAR. 

We further verify that the energy dependence 
of the indirect electron longitudinal shower infor- 
mation behaves as expected from electromagnetic 
shower models. 

2. The BABAR Detector 

The BABAR detector consists of a tracking sys- 
tem surrounded by a dedicated P1D device, a crystal 
calorimeter, and an array of flux return iron plates 
instrumented with muon detectors. A detailed de- 
scription of the BABAR detector can be found in 
Ref. [2] , here we briefly describe those systems im- 
portant for the algorithms described in this paper. 

The BABAR tracking system is composed of a sil- 
icon vertex tracker (SVT) comprising five layers of 
double-sided detectors and a 40-layer central drift 
chamber (DCH). Operating in a 1.5 T solenoidal 
magnetic field, it provides a transverse momentum 
(Pt) resolution 



a PT 

p T 



(0.13 ± 0.01)% • p T © (0.45 ± 0.03)% (1) 



for detecting charged particles, where pt and its 
rms error ap T are measured in GeV/c. Both the 
DCH and the SVT measure the specific ionization 
(dE/dx) of charged particles which pass through 
them. The dE/dx resolution from the DCH varies 
as a function of track polar angle, having an aver- 
age value of around 11 % for the majority of physics 
processes that we study in BABAR. The dE/dx reso- 
lution from the SVT is typically about 16 %. 

The BABAR ring-imaging Cherenkov detector 
(DIRC) provides dedicated charged particle iden- 
tification in the central part of the detector. The 
polar angle coverage in the laboratory frame is 



—0.84 < cos 8 < 0.90. The Cherenkov angle resolu- 
tion of the DIRC is measured to be 2.4 mrad, for 
the quartz refractive index of 1.473, which provides 
better than 3 a separation between charged kaons 
and pions over a broad kinematic range. 

The BABAR electromagnetic calorimeter (EMC) 
consists of an array of 6580 CsI(Tl) crystals, which 
encloses the tracking system and DIRC. The crys- 
tals have a truncated trapezoidal shape, and are 
finely segmented in the plane transverse to parti- 
cles coming from the nominal e + e~ interaction point 
(IP), with a typical cross-section of 4.7 x 4.7cm 2 at 
the front and 6.0 x 6.0 cm 2 at the back. The crys- 
tals range in depth between 16 and 17.5 radiation 
lengths (the radiation length of CsI(Tl) is 1.85 cm), 
with the crystal axis pointing back roughly to the 
IP. The EMC geometry can be approximately de- 
scribed as a central cylindrical barrel, divided into 
forward (26.93° < 8 < 90°) and backward (90° < 
8 < 140.81°) regions, plus a conical forward endcap 
(15.76° < 8 < 26.81°). The crystals are staggered 
so that their front face presents a nearly normal sur- 
face to particles coming from the IP. The EMC cov- 
ers about 90 % of the polar angle and all of the az- 
imuth in the center-of-mass system of the collisions 
produced in PEP-II (see Fig. 1). 



229 5 External 
Support 




Interaction Point 



Fig. 1. Longitudinal cross section of the EMC showing the 
top half of the detector. All dimensions are given in mm. 

The energy resolution of the calorimeter has been 
measured in-situ using a combination of radioactive 
sources, symmetric decays of tt° and r\, and Bhabha 
events, and can be described as: 



<j e _ (2.30 ± 0.30)% 
~E~ ~ 



VE(GeV) 



(1.35 ±0.22)%, 



(2) 



where E and <je refer to the shower energy and its 
rms error, measured in GeV. The angular resolution 
is limited by the transverse crystal size and the dis- 
tance from the IP. It can also be parameterized as 
an energy-dependent function 
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The EMC is surrounded by a series of iron plates 
arranged as coaxial octagonal cylinders about the 
BABAR symmetry axis. These plates form a high- 
susceptibility path for the magnetic flux generated 
by the BABAR solenoid to close on itself. Between the 
iron plates are resistive plate chambers and limited 
streamer tubes with binary readout, used to track 
muons and provide crude neutral hadron detection. 
The innermost layers of muon chambers act effec- 
tively as a 'tail catcher' for the EMC, detecting par- 
ticles from showers that leak out the back. 

3. Longitudinal Shower Depth Variable 

To derive longitudinal shower development infor- 
mation from BABAR we exploit the fact that most 
particles do not enter the calorimeter exactly par- 
allel to the crystal axes. A non-zero entrance an- 
gle transforms the transverse crystal segmentation 
into an effective longitudinal segmentation, provid- 
ing some depth information. Because the effective 
longitudinal segmentation is poor (often fractional) 
and different for every particle, we do not attempt a 
full parameterization of the longitudinal shower de- 
velopment. Instead, we characterize the shower by 
the first moment of its longitudinal development, 
which we call the Longitudinal Shower Depth (AL). 
The AL value is closely related to, but not identical 
to, the position of the electromagnetic shower max- 
imum, as is discussed in Appendix A. 

The AL variable is a geometric quantity which 
exploits the fact that the track and the cluster both 
sample different two-dimensional projections of 
the three-dimensional shower spatial distribution. 
When the track direction is not parallel to the crys- 
tal axis, these projections are not fully degenerate, 
and they can be combined to extract the otherwise 
unobservable, third (longitudinal) dimension. 

Three effects are responsible for the fact that the 
track direction and the crystal axis are not collinear. 
First, the magnetic field bends the track as it passes 
through the tracking volume. Second, the width of 
the beamspot in the beam direction causes tracks 
from the IP to have a different polar angle from 
that of the axis of the crystal they strike. Finally, by 
design, the crystal axes of the BABAR calorimeter do 
not project perfectly back to the nominal IP, which 
reduces the chance of particles from the IP passing 
perfectly between crystals. 



As part of computing AL we describe the 
calorimeter cluster as a directed line segment in 
space. We first compute the two-dimensional clus- 
ter centroid using the standard BABAR algorithm, 
which takes the weighted average of the crystal 
center positions at a nominal depth of 12 cm, using 
a logarithm of the crystal energy as weight [4]. We 
then compute the weighted average direction of the 
crystal axes, using the energy in each crystal as 
(linear) weight. The cluster line segment is defined 
to pass through the cluster centroid, and point in 
the average crystal direction. The starting point 
of the cluster line segment is taken as the average 
position of the crystal front faces projected along 
the average direction. 
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Fig. 2. (color online) Schematic view of how AL is calculated. 

We then calculate the point of closest approach 
(POCA) in three dimensions between the extrapo- 
lated track trajectory and the cluster line segment, 
using an iterative algorithm. The POCA is the point 
where the track and cluster projections of the par- 
ticle trajectory are most consistent. 

Conceptually, we define AL as the path distance 
the track travels in the calorimeter's sensitive vol- 
ume in reaching the POCA. In practice, we define 
AL as the distance along the cluster line segment 
of the POCA, divided by the cosine of the angle be- 
tween the track direction and the cluster line seg- 
ment direction, given algebraically as: 



AL 



(^POCA — ^Front) \ ^Cluster 
COS/3 



(4) 



where rpocA is the position of the POCA, rp ron t is 
a point on the front face of the crystal, fciustor is a 
unit vector in the direction of the cluster line seg- 
ment, and (3 is the angle between the track direction 
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and the cluster axis direction. This quantity approx- 
imates the material path distance, but is much sim- 
pler to compute. Our definition of AL ignores the ef- 
fects of track curvature and crystal-face staggering, 
which are negligible on the scale of the resolution 
we achieve on AL. The definition of AL is presented 
graphically in Fig. 2. 

4. Particle ID Control Samples 

We evaluate the usefulness of AL for particle iden- 
tification with the control samples of relatively pure 
electrons, pions, kaons, muons and protons. The con- 
sidered control samples are selected from the data 
collected with the BABAR detector at the PEP-II 
asymmetric energy (3.1 GeV on 9.0 GeV) e + e~ col- 
lider, operating near the T(4S) resonance which sub- 
sequently decays into a BB meson pair. In the fol- 
lowing, we briefly outline the salient features of the 
control samples. 

We select electrons^ in both radiative and non- 



radiative Bhabha events [e + e 



"(7)] by 



utilizing requirements on the energy deposit and 
shower-shape variables in the EMC, and by reject- 
ing track candidates consistent with being muons. 
Based on Monte Carlo studies, the purity of this 
sample is found to be 99.9 %. As Bhabha events 
provide mostly high momenta electrons, we con- 
sider the two-photon mediated process e + e~ — > 
(e + e~)7*7* — > (e + e~)e + e~ to enhance statis- 
tics in the low-to-medium momentum range (p < 
3 GeV/c). The selection requirement for this pro- 
cess provides a clean sample of electrons with purity 
comparable to that of the Bhabha events. 

Pion candidates are selected from the decay pro- 
cess K® — > 7r + 7r~ and the e + e~ — > t + t~ events 
with 3-1 track topology. The purity of the sam- 
ple selection, determined with a mixture of simu- 
lated BB decays and e + e~ — > qq continuum events, 
is found to be ~ 99.5 %. Pions from r-pair events 
are affected by kaon contamination, having a purity 
of about 97%. However, this sample provides high- 
momentum pions not found in the Kg — > ir + ir~ sam- 
ple. The electron contamination in the r-pair pion 
sample is very small. 

In addition to the electron and pion control sam- 
ples, we use kaons selected from the decay chain 
D*+ -> D Q Tr+,D a -> R-T7+; protons from the A 
decay, A — ► p7r~; and muons selected in radiative 



muon-pair events, e + e~ — > p + /i — 7. These samples 
have the best purity for the corresponding track- 
candidate selection in BABAR, which is comparable 
to that of the considered electron and pion samples. 

5. Electron ID Performance 

To test the impact of AL on electron identifica- 
tion, we start with an electron selector based on two 
standard variables: the ratio of the shower energy 
deposited in the calorimeter to the momentum of 
the track associated with the shower (E/p) and the 
lateral shower moment, defined as 



LAT = 
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Here N is the total number of crystals associated 
to a shower, Ei is the energy deposited in the z-th 
crystal such that E\ > E2 > ■■ > En, ri the lat- 
eral distance between center of the shower and i- 
th crystal as defined earlier, and tq = 5 cm which 
is approximately the average distance between two 
crystals. Like AL, these variables involve calorimc- 
try and tracking measurements only. By evaluating 
the incremental improvement given by adding AL 
to this selector, we test the impact of AL including 
the effect of possible correlations between AL and 
other similar variables. 

We describe the electron identification perfor- 
mance in terms of the electron efficiency and the 
pion misidentification of the algorithm. We use 
the TMVA (Toolkit for Parallel Multivariate Data 
Analysis) package [5] to build a global likelihood 
function using E/p and LAT, together (or not) 
with AL. 

To study the pion misidentification probability for 
a given value of electron ID efficiency, first we define 
a likelihood ratio, Rl, for each track candidate in 
the considered signal and background samples by: 



Rl 



L.c 



(6) 



Here, the signal and background likelihoods (Lg, 
Lb) are products of corresponding probability den- 
sity functions (ps, Pb) of the three discriminating 
variables: 



Electrons denote both electrons and positrons. 



Ls(i) 



Y[ps,j(i)- 

j'=i 



(7) 
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After that, for a given value of the likelihood ra- 
tio, the signal identification efficiency and the back- 
ground misidentification probability are calculated. 
This is done for different momentum bins, separately 
in the forward and the backward barrel, and the end- 
cap regions, using the control samples of electrons 
and pions discussed in the previous section. 
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Fig. 3. (color online) Pion misidentification probability 
as a function of electron identification efficiency in the 
forward Barrel region for a typical low momentum bin: 
0.2 < p < 0.4GeV/c. 

Table 1 

Comparison of pion misidentification probabilities at 90 % 
electron identification efficiency in the case where the likeli- 
hood function is defined with (or without) AL. 
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Figure 3 shows the electron efficiency vs. pion 
misidentification probability for a typical low mo- 
mentum bin (0.2 < p < 0.4GeV/c) in the forward 
barrel EMC. It is evident that for any given value of 
electron identification efficiency the likelihood func- 
tion based on AL gives a lower pion misidentifica- 
tion compared to the case where AL is not included. 
Table 1 summarizes results obtained across the full 
kinematic range for various parts of the EMC. There 
is a clear improvement in the performance for the 



backward and forward barrel regions, while for the 
endcap region (where high momenta particles are 
mostly abundant), we find a marginal improvement. 
This is because the discrimination power of AL di- 
minishes with increasing energy. 

6. Charged Particle ID Performance 

The AL variable can also be used to enhance gen- 
eral charged particle identification, as it is sensitive 
to the differing longitudinal shower development of 
different particle types. This is demonstrated in Fig. 
4, which plots AL for different species of particles, 
broken down into four track-momentum bins, using 
the BABAR data PID control samples described in 
Section 4. 

Figure 4 shows a clear distinction between the AL 
distributions of different particle species, particu- 
larly for momenta below 600 MeV/c. These distribu- 
tions can be basically understood in terms of the dif- 
ferent energy loss mechanisms at work; for instance, 
low momentum protons are highly-ionizing, and so 
deposit most of their energy early in the crystals. By 
contrast, electrons deposit their energy near shower 
maximum (roughly 10 cm), while muons with mo- 
menta above 200 MeV/c are minimum ionizing and 
so distribute their energy uniformly along their path 
through the EMC. Finally, pions and kaons produce 
broad AL distributions, corresponding to the large 
variability of hadronic showers. 

In this momentum range AL provides 0.8 a pion- 
muon separatiorPI. compared to 1.5 a separation 
from the DIRC, less than 0.1 a separation from ei- 
ther DCH or SVT dE /Ax [8] , and essentially no sep- 
aration from E/p. Thus AL provides an useful cross- 
check to the DIRC when identifying muons at these 
momenta, and provides the best muon-pion separa- 
tion for the 15 % of the BABAR solid angle covered 
by tracking and calorimetry but not by the DIRC. 

At momenta above 1.25 GeV/c, the decrease in 
magnetic bending reduces the angle between the 
track direction and the crystal axis, degrading the 
resolution of AL. Additionally, the longitudinal pro- 
file of energy deposition for different particle types 
tends to converge in this momentum region. Some 
separation power still comes from different widths 
of AL distributions for electrons compared to other 
particles, but this is a weak discriminant compared 

2 We define a as the difference between the average of the 
muon and pion AL dis tributions div ided by the quadratic 
average of their rms, yj (cr 2 , + <r 2 )/2. 
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Fig. 4. (color online) Distributions of AL for different types of particles in different momentum bins, selected from the BaBar 
data control samples. Note the differences in the x-axis range. Each histogram has been normalized to unit area, to better 
show the AL distribution shapes. 



to other PID variables available in this momentum 
region. 

The impact of AL on muon identification at 
BABAR has been evaluated using a powerful muon 
selection algorithm which uses a bagged decision 
tree [9] to combine many input variable^]. This 
algorithm was trained and evaluated using indepen- 
dent subsets of the data control samples described 
in Section 4. Compared to an older algorithm which 
does not use AL, the minimum muon momentum 
for which the selector has at least 50 % efficiency 
(at a fixed pion misidentification probability) was 
reduced from 800 MeV/c to 270 MeV/c [8]. This 
improvement in low momentum muon selection ef- 



3 This muon selector has 30 input variables. 



ficiency is expected to have a significant impact on 
several important BABAR physics measurements. 

7. Energy Evolution of Shower Depth 

Due to the exponential nature of electromagnetic 
shower development, we expect the shower maxi- 
mum position to evolve logarithmically with energy. 
A standard model describing the logarithmic energy 
dependence of longitudinal development is given in 
Ref. [6]. Because AL is closely related to shower 
maximum, we also expect that to evolve logarithmi- 
cally with energy as well. 

Figure 5 shows AL for electrons against the nat- 
ural log of the electron energy in the BABAR data 
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Fig. 5. (color online) a) Dependence of AL on energy for 
electrons in the BaBar. data two-photon and Bhabha control 
samples. The points represent the mean from a Gaussian fit 
to x-slices of the full histogram. In b) The same points are 
shown with an expanded y scale. The overlaid fit is to a 
shower model function described in the text. 

two-photon and Bhabha control samples described 
in Section 4. The electron energy is estimated as 
the track momentum calculated at the point where 
it enters the EMC. The two-photon data cluster at 
low energy, the Bhabha at high. The points are from 
slicing the data along the horizontal axis, fitting 
each slice to a Gaussian, and plotting the Gaussian 
mean value with its error. The points are fit to a 
function based on the shower model parameteriza- 
tion of AL for electrons described in Appendix A. 
This function has two free parameters: the material- 
dependent shower scale factor plus an overall 
shift in AL. The fit quality is reasonable, showing 
consistency of our data with the model. The shift 
value we obtain of 0.3 cm is consistent with that ex- 
pected due to the roughly 0.2 % per cm light atten- 
uation measured in the BABAR CsI(Tl) crystals. The 
value of the shower scale factor we obtain is b = 
0.492 ± 0.001, in good agreement with the value ex- 
pected for CsI(Tl) of 0.492 ± 0.010, where the error 
on the expected value comes from interpolation in 



atomic number of the data provided in Ref. [7] , plus 
the small residual energy dependence of this param- 
eter over the relevant energy range. 

The points in Fig. 5 show both the statistical 
(small bars) and total (statistical plus systematic, 
larger bars) errors. The dominant systematic er- 
ror comes from the alignment of the EMC with re- 
spect to the BaBar tracking system, which we es- 
timate has an uncertainty of 0.2mrad in the rota- 
tion around the beam axis. Because this misalign- 
ment can depend on the track polar angle, and be- 
cause the boost of PEP-II correlates energy with 
polar angle, we include this systematic as a point- 
by-point error. Other systematic effects are from 
hadron backgrounds in the sample, the uncertainty 
in the amount of material before the calorimeter, 
and changes in light attenuation in the crystals due 
to radiation damage. These latter are largely accom- 
modated by allowing the depth offset parameter in 
the fit to float. 

8. Monte Carlo Performance of AL 

We have also studied AL in Monte Carlo events 
simulated in the BABAR detector, and compare it 
with real data. This provides a useful test of the 
shower development model used in simulation. 




AL (cm) 

Fig. 6. (color online) Distributions of AL for different 
types of particles in the intermediate momentum bin 
0.4 < p < 0.6GeV/c produced in Monte Carlo BB events 
processed by the BaBar. detector simulation. The histograms 
are normalized to unit area to better display their shapes. 

In Fig. 6 we show AL distributions for differ- 
ent species of particles selected in simulated control 
samples for the intermediate momentum bin (0.4 < 
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p < 0.6GeV/c), chosen because it shows maximum 
separation among various particles as discussed in 
Section 6. The Monte Carlo AL shapes are very sim- 
ilar to those observed in BaBAR data (see Fig. 4) , in 
this and other momentum regions. 

Furthermore, we have verified the logarithmic de- 
pendence of AL on the incident electron's energy 
in simulated samples of electrons from Bhabha and 
two-photon processes. Fitting the AL vs. ln(L) dis- 
tribution of these simulated electrons to the shower 
model function described in Appendix A gives a 
value for the shower scale factor b — 0.493 ± 0.001, 
consistent with expectations for a CsI(Tl) calorime- 
ter and with our measurements in BABAR data. 

9. Conclusions 

In this paper, we have presented a technique for 
extracting longitudinal shower development infor- 
mation from a crystal calorimeter in conjunction 
with a precision tracking system. When the derived 
quantity AL is used in electron identification, we 
have shown that the algorithm performance is sig- 
nificantly enhanced, especially in the low momen- 
tum region. We have also shown that AL can signifi- 
cantly improve particle separation for other types of 
particles, particularly muons and pions, in the low 
momentum region. By studying the energy evolution 
of AL for electrons in BABAR, we have established 
that it behaves consistently with expectations from 
a standard longitudinal shower model computation. 
We find that the BABAR simulation well reproduces 
the observed AL behavior. 
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Appendix A. Energy-weighted shower depth 

The mean longitudinal profile of the energy de- 
position in an electromagnetic shower can be well 
described by a gamma distribution [6] : 

^§=E f(t) with/(t) = f^y (&*) a-1 be~ bt ,(A.l) 

where t is the distance measured in units of radiation 
length Xo, Eq is the energy of the incident particle 
and E is the energy deposited by the particle at a 
certain distance t. Here b is a material-dependent 
shower scale factor [7], and a = (1 + bt max ), where 
^max denotes the shower maximum and is expressed 
as follows: 

t mm = hi(Eo/E c )+C j , j = e,7, (A.2) 

where C e = —0.5 for electron-induced showers and 
C 7 = +0.5 for photon showers. To a good approx- 
imation, the critical energy E c is given as E c = 
0.8 GeV/(Z + 1.2) = 0.0145 GeV [for CsI(Tl), aver- 
age Z = 54]. 

We can then solve for the energy-weighted lon- 
gitudinal shower depth, AL, by integrating the 
fractional energy deposition per radiation length in 
Eq. (A.l), f(t), with the corresponding path length 
t: 

L L 

AL = jt f(t) dt = Y^)j ( bt ) a dt ( A - 3 ) 
o o 

Here L indicates the average length of the CsI(Tl) 
crystal in units of Xq and varies between 16 to 17.5 
for various regions of barrel and endcap EMC. For 
simplicity, we took L to be 17. 

We note that the AL expression in Eq. (A. 3) has 
a roughly logarithmic energy dependence [see the 
definition of a and Eq. (A.2)], and that the shower 
scale factor b is roughly independent of energy. We 
can therefore approximate the energy and material 
dependence of AL by Taylor expanding Eq. (A. 3) 
in terms of ln(E) and b. Keeping only the leading 
linear terms: 

AL « 10.4 - 4.9A6 + (0.9 - 66.0A6) ln(L), (A.4) 

where Ab is the deviation of b from a nominal value 
of 0.5. This expression can be used to describe the 
AL dependence on ln(L), for material with effective 
Z values near that of CsI(Tl). 
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